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ABSTRACT 


Computed Channel 2 or window radiances for the Nimbus II 
optical-sensor system were kindly made available to the author by 
the Meteorological Satellite System of ESSA for each of the 106 
model atmospheres considered by Wark et al. (1963) 

Multiple regression equations relating the emitted black-body 
intensity from "black" interfaces to the independent variables (1) 
"observed" specific filtered radiance and (2) total optical path, 
were set up for each of seven zenith angles in the range © =0° to 
60°, and for each of 62 randomly selected cases from the Wark sounding- 
catalog. Both variables gave high statistical significance, with the 
former accounting for the primary part of the variance of the dependent 
variable, but with variable (2) always contributing in such a way as 
to account for some of the atmospheric absorption along the sensing 


direction. 
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1. Introduction 

Wark et al. [7a] on and other authors have shown that both the 
Surface air temperature under clear sky conditions, and the cloud-top 
temperature when the sky is covered by relatively dense clouds, as 
Statocumulus or altostratus, can be estimated using outgoing terres- 
trial radiation intensity observations obtained by Tiros II, III and 
IV channel 2 radiometers. Channel 2 radiometers are filtered to 
receive only infrared radiation in the "window" region of the 
Spectrum, and in the Nimbus II this window region was 10-1144 
Usually the estimated temperature was several degrees too low under 
clear sky conditions. 

By utilizing Nimbus II channel 2 radiometer computations of 
sensed outgoing terrestrial radiation intensity (as computed by the 
Meteorological Satellite Laboratory) for 62 atmospheric models -- 
all of the 47 cloud-undercast cases, and 15 randomly selected clear- 
sky cases, multiple linear regression equations are developed to 
solve for the angular dependence of total outgoing specific intensity 
1(@ ) values. The multiple linear regression analysis uses the 
BMDO3R program as adapted to the CDC 1604 computer. 

In order to describe radiative temperatures at the earth's 
Surface under clear sky conditions, and at the cloud-top for under- 
cast conditions an estimate of the total specific intensity as viewed 
by Nimbus II must be made as a function of the emitted black body 
intensity and of an absorptivity dependent upon the effective optical 
path (LASEC QO ). Since, in general, window radiation is more 
readily available, it was decided to use specific filtered (window) 


radiation, rather than total specific intensity. The latter may then 


be inferred as a function of interface temperature and optical path, 
as is done in Section /. 

Graphs/tables are then constructed which allow simple compu- 
tation of the surface air temperature under clear skies, and cloud 
top temperatures, using total water vapor mass, zenith angle, and 


satellite radiometer measurements in Channel 2. 


2. Data Sources 

The data utilized consisted of the outgoing terrestrial radiation 
intensity computed by the Meteorological Satellite Laboratory from 
window region(Channel 2) data obtained from the Nimbus II satellite. 
106 atmospheric models were developed by the Meteorological Satellite 
Laboratory in order to calculate the outgoing intensities rae of 
which 62 selected cases are used in this paper. The 62 cases are 
composed of 4/7 cases undercast at various altitudes and 15 clear 
cases, which are identified in Appendix I. The models developed 
by the Meteorological Satellite Laboratory were taken from radio- 
sonde data over all latitudes and seasons. All the soundings went 
to at least 25 mb, and the temperatures were extrapolated parallel 
to a standard atmosphere to the pressure of 0.Ilmb. Stratospheric 
hunidity was extrapolated to conform in general to standard-atmosphere 
values. The sky condition, the temperature at selected levels and 
the temperature of the undercast level, as well as the total water 
vapor and oZone masses were also listed as part of the sounding 
catalog in ial) 

In order to deduce total emergent flux and the angular dependence 


of these flux estimates, the specific filtered intensity in the 10-1] 


band for each of the 62 atmospheric models was used at each of the 
four given zenith angles ( 6 = Oo: 20°, 459 and 60°). The inten- 
sities at three additional zenith angles (© = 109, 32.5° and 52.59) 
were obtained by interpolation. 

The radiation intensities of the 62 models were given in ERG/ 
(cm2 SEC STRDN) and are converted to watts per meter“ per ‘strdn 


in this thesis. 


3. Nature of the Data, and of the Model 
The Nimbus II satellite uses both a High Resolution Infrared 
Radiometer (HRIR) scanning the 4.0 micron window, and a Medium 
Infrared Resolution Radiometer (MRIR) which scans the infrared 
Spectrum in regions where the far infrared region of the solar 
Spectrum gives only insignificant contributions. Only Channel 
2 radiation of the MRIR segment is used, and the sensed radiation 
comes from an area of approximately 40 miles across in the hori- 
zontal (based upon a half-cone viewing angle of 2.9°). The solid 
angle thus presented at the sensor is sufficiently small so that each 
measurement must be considered a specific intensity rather than flux. 
The emergent flux is a more important quantity than that of 
outgoing specific intensity from the meteorological point of view, 
because the latter applies to a specific direction O only, whereas 
the emergent flux gives total radiative heat loss from a column. 
In atmospheric radiation problems it is convenient to treat the 
ground and cloud surfaces as black bodies | 11 . As shown in 


L3, pp. 81-82 | the intensity and flux of radiation from a black body 


are respectively 


¥ 
T,-—t AND Fp =OT (1) 


in terms of the black-body surface temperature T,. In (1), the 


constant UO of the Stefan-Boltzman law has the value 
~7 P 
a =0.5664 X/0 WATTS /m 


so that all values of FR are in watts /m* and those of Ip are in 
watts /m2/steradian. 

The emergent filtered intensity, I © (©), scans primarily 
the 10-1144 window region. The nature of the filter-transmissivity 
function in this band is an important consideration. The techniques 
by which the filter-functions are employed in conjunction with the 
sensor absorptivity is discussed in considerable detail by Wark 
et al i p. 3, and Supplement, p. Dale These same considerations 
have been applied to the measured values of l@ (©) for Nimbus II, 
and the computations provided by the Meteorological Satellite Lab- 
oratory are presumably available to other interested workers. The 
interesting facts are that (1) the thermister sensitivity is not 
unity in any of the wave band intervals, (2) that window filter 
function is not unity over the window, nor is it zero outside this 
spectral region. However, suffice it to say Wark's filtered radiance 
(7, Eq. 3 | multiplied by the effective thermistor emissivity 
provides the "observed window radiation" of this study. Such data 


is of the form 


—_—_—_// 


/0 
eS dT, 
Z, (6) C( He (2) 


with C known as a quasi-constant with respect to l@ (@) over a 
wide angular range. Theoretically one could then enter a table of 
I @(@) FREES = ONY IGE by C, and recover the appropriate temper- 
ature for (2a) from the Planckian monochromatic intensity function. 
However, even in the so called window region there are small atmospheric 
depletions due to both water vapor and ozone, which should be restored 
if accurate "black-surface" temperatures are to be estimated. Thus, 
in (3) below a depletive effect is added back to recover part of 
the initial beam intensity. 

Hence the linear regression was set up with a heuristically 


physical basis, as follows: 
2 
if 
=e = Ar A,Lo(O)lI+Az 4(=) SEC 6 | (3) 


Here * has been multiplied by a single temperature-ratio correction 
factor of the Lorentz-type, since the individual computations of U, 
did not include such factors in the layer-by-layer summations. The 
inclusion of this temperature-ratio term, along with the constant Ag, 
(to be determined) was aimed at inducing a more homogeneous strati- 
fication of the extinction treatment, since the cloud top cases 
investigated here varied considerably with respect to pressure level. 


Thus, (3) becomes 
Y Te I 
— =A, tA,Lo(6)+42A (Z) LQ) ysece (4) 


Moreover in view of the one-to-one (expected) relationship between 
l ® (CE) and Gaal" remarked upon above, a linear statistical 


relationship of form 


1] 


- 1) (5) 
I, (0) = deg ( GT tiles 
is proposed, with both ky and ko quite small. In fact, apart from 
sensor-noise, and second-order depletive effects, ko should be 
negligibly small and is treated thus here. (4) then assumes the 


form 
@ 
TTA TA, (nIo( Or TALA, hk) r*(Z) SECO (6) 


where Aj, Ay; rime and k would each have to be determined by standard 
best-fit methods. The grouping nm o ( e ) is emphasized in order 
to show the conversion of intensity to flux units throughout (6). 


We therefore arrive at an alternative form of (6) 


Y, = 5+ B, X2+ & Xs = 
where the independent variables are 
X2 = 71, (0) vs 
Kg =o79(Z)*y, SECO 
and the dependent variable is 
Y,=oT" (8b) 


The aim of the working equation, (7), is to determine how well 


the effective black-body temperature can be described statistically 


by use of observed values of filtered radiance l@ ( =) and of a 
more complex variable Xa which appears to depend upon the dependent 

, 2 
variable, but actually its variance, Ts a= CONST (Ty sece). 
However X3 was formed as shown with values of UL also Known in advance 


for each of the 62 soundings used, in accordance with 


Rureace 
: 
Ue= “5 J t (Gas) §P ’ 
P=0./ug 

The choice of the standard temperature le = 300°K was used in 
connection with the full Lorentz-type broadening factor as set back 
into (3). This value was adopted based primarily upon experimental 
transmissivities of Palmer [5] who found that he could reproduce 
the theoretical absorption coefficients of Yamamoto as listed inf 7], 
computed for T = 260°K. This difference is well brought out by Wark 
et al. [ 7, eq. 5 | , which suggests that in order to retain U, in the 
general format of (9), the exponent of (p/p) should be slightly 
less than unity. But then if strict conformity with Lorentz-type 
broadening is required, the standard temperature IT, is raised to 


300K. 


4. Data Processing 

As data samples for the testing of (7) at seven different angles 
in the range O= 0° to 60°, the cloud-top and/or surface temperatures 
Wegeslce@mds INputsfor 1%) = iy for each of the 62 soundings chosen 
from M.S.L. Report Number 10 [7] . For each such sounding, values 
of I © (0) at angles S) = Oye 209, 45°, 60°, were also given for 


Nimbus II. The I (0) sample was arbitrarily enlarged by linear 
interpolation to corresponding values at iy 32.59, 52, 5- giving 
seven angles for each sounding with seven corresponding values I@ (O). 


The variable 


= 1E 
Xe =(TT") (sb) ly SECO 


was formed by transgeneration using a program available for this 
purpose at U.S. Naval Postgraduate Schoo] Computer Facility. There 
were 7 collections of punched data cards for each sample set (Y),%X%o,X9)5 
each collection corresponding to a different zenith angle applied to 
the set of 62 soundings. 
The multiple regression program BMDO3R, as adapted to the 
CDC 1604, was utilized to solve for the value of the best fit 
coefficients for each of the seven zenith angles. The total input 
to the program consisted of /x62 or 434 data cards. Table 1 shows 


a sample input for one zenith angle. 


5. Statistical Results and Inferences 

The multiple linear regression equations obtained from the program 
are shown in Table 2, in the form of (7) and (8). It can be readily 
seen by examining the formulas of Table 2, that the regression coeffi- 
cient of the first independent variable considerably over-weights 
that variable when compared to the regression coefficient of the 
second variable, the depletion factor. 

In order to examine the question of a significant specification 


of oT? by the independent variables X. and X3= oT "(4)" U, secO , 
0 


Table 1 


Sample-set Employed in the | 
Multiple Regression Analysis * 


2) Xo = lo (0) X,(watts/Mé) 4 (watts /M*) 
10 D7 Ie 40.084 262.854 
10 22.008 108.177 323). 905 
10 07a 16.523 2/5. beg 
10 26.423 495.751 373.705 
10 paelsy, 11.423 220.381 
10 sey on 40.668 247.104 
10 27.962 1071.279 406.182 
10 16.987 50.299 262.854 
10 24.897 Sieh || 16) SOG 70C7 
10 21.474 129.534 519209 
10 a pg Wake 181.637 563707 
10 opal) Sis 1.694 158.511 


*The above 12 samples represent 12 of 62 samples for the zenith 


angle © = 10°. 


Es 


0° 
10° 
20° 
32.50 
450 
52.59 
60 


Table 2 


Multiple Linear Regression Equations Obtained from 


the BMDO3R Program for Seven Zenith Angles 


2602/0 
58.869 
Sama 2 
629260 
56.224 
Symes 
63,687 


. 3788X 


. 1866X 


.9226X. 
.8840X., 


.9233X, 
.2479X> 


.1443X, 


.00724X 


.00983X 


.01121X 


.01001X. 
.01084X, 


3 


.00998X. 


c: 


.01125X3 


3 


the F=statistics in Table 3 will be used as obtained from the output 
of the multiple linear regression program. 

The BIMDO3R program tests for significance of variables 
introduced into a multivariate linear regression equation. The most 
readily adaptable test, in the analysis of variance, is the so-called 


F-test where the 


F-estatistic = Mean squares explained by predictor (10) 
Mean squares unexplained 


In this experiment Xo is the primary specifying variable, and for 


its test F = FS 60 is defined as 


F - Sum of squares explained 
1,60 [Sum of squares unexplained | 
(N-2) 


The subscript numbering (1,60) denoted the number of degrees of 
freedom associated with numerator and denominator, respectively. 

The addition of the depletion factor, symbolized by X3, requires 
the definition of a new F-statistic specifically restricted to the 
testing of X3, after the contribution of Xo has been excluded. This 


statistic will be denoted F\ and be defined as 


- Additional Sum of squares explained any) 
: 59 pene of squares unexplained 
The critical value at the 99% confidence level of belief is 7.08 for 
both F and Sa and these critical values are surpassed at the 99% 
confidence level for each variable added, with the exception of 


X3 at ® = 20° and B2e5 However the ratio 


Teo) = o(10%) & 0 (10%) 
j 


7,0 


TABLE 3 


Some Statistics Relating to the Multiple Linear 
Regression Equations for Seven Zenith Angles 


Partial Mult Computed 
QO vari- Mean Reg.Coeff Corr Coeff Corr Coeff F & F' Values 


able 
0° =X, ~——*19..690 11.92261 .99669 23927 .718 
Xz 272.446 .01001 = . 55242 25.912 
Y; 296.471 9987 
10° Xy_—s«d19. 681 11.88413 —.99463 14765 .000 
Xz 275.661 .01084 .49824 18.482 
Y; 296.471 .9980 
208 Xo leee2e 12.37883  .96464 2088.534 
X§ 289.959 .00724 =. 14609 1.287 
¥,  290ma7) 9859 
32.59 X, 19.436 11.92329  .96001 2972.724 
X3 323.062 .00998 24503 3.769 
Y, 296.471 .9876 
45° =X, 19.248 12.24789 .99519 22457562 
Xz 385.321 .00983. = .61318 35.550 
Y, 296.471 .9982 
52.59 Xo ~=—«19.131 12.18661 .98884 7231.164 
X3 458.220 01125 =. 56546 27.733 
Y) 296.471 .9959 
609 xX,  =+18.720 12.14434  .94758 1427.700 
Xz 544.308 .01121 = .36280 8.943 
Y, 296.471 .9798 
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and is so far beyond critical that the contribution of Xo is highly 
Significant, and certainly cannot be accounted for by a chance 


contribution in any of the seven cases. The contribution of Xo 


as test by F 1-59 (Xs) 1 
.08 
in all cases except @ = 20° and 32.5°. 

It should be noted that the coefficient of the depletion 
factor X, is positive in all cases; that is, it helps to account 
for a small but significant part of absorption from the originating 
flax Cai. which part has been absorbed in the atmosphere. It has 
therefore been tentatively concluded that all 7 angular regression 
formulas afford significant information with respect to the atmo- 
spheric depletion, as well as to the primary variable 1@ (@). 
Here, in effect some weight is being given to the sign test of the 
regression coefficient of Xo even for © = 20° and 32.5°, where the 
F' test falls short. The F-test failure in the case @ = 32.5 is due 
to the inclusion of O = 20° data in the formation of the O = 32.5 
Sample. 

An explanation for the lack of significance of the depletion 
factor for 6 = 20° could lie in the interplay of water vapor and 
ozone radiations, weak as both of these are likely to be. 

The shortcoming of the statistical model must be recognized: 
only one absorbing element has been introduced, namely water vapor, 
which absorbs weakly in the "window" region| 7, mig. 3 | of the 
atmosphere. In addition, the small ozone content introduces a 
moderate absorbing band where the water vapor concentration is 


sparse (15 km and above). The emergent intensity formula for a 
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Single-absorber atmosphere is given by 


z, (6)=Ley (Trop) 2 ory ( UgSECO) 928/(ar) (12) 


and, in a wave number interval where two absorbers give overlapping, 
the transmissivity ty is generally taken as the product of the 
individual transmission functions for the gases taken separately. 

First of all it appears that in the direct computation of top- 
transmitted radiance, there may be an angle where overlap-effects 
give larger discrepancies than at extremely close or extremely 
distant angular ranges. Secondly, the statistical attempt of 
synthesizing two absorbers by putting in just water vapor could 
not account for all of the variance attributable to depletion. 
Minor variations in the proportionality constant C of (2) undoubtedly 
do occur for any given setting of filtering wave band limits. 

On the whole, however, a physically-based statistical mode] 
which has a multiple correlation coefficient 4 O98 1S qutee 
adequate for specification of the dependent variable, and in the 
present case, of accounting for approximately one-third of the 
unexplained variance left over when I (6) alone 1S used in 


specification. 


6. Estimation of Surface and Cloud Top Temperatures 

Utilizing the multiple regression equations listed in Table 
2, for the seven zenith angles, Fig. |] and Appendix II have been 
developed to afford close approximations of the surface and/or 


cloud-top temperatures. The basic input data for use in these 
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graphical aids are the Nimbus II Channel 2 radiation measurements, 
together with the water-vapor path in direction 6 . These graphs 
are based upon (15) which is derived below. 


Since by (7) 


TT1= B + &XxgrBghiz 


we may define, a unique black-body flux using the mean values, 
B, and Bo: 


TT! = +B Xz 


7 = ica ] 


“™ 
Then making the approximation? T° = ope in X3, we have 


a if 
Ay “A R (14) 
gTt=07, + 8|(o73)(Ze) 4g Sec | 


By use of the binomial expansion carried out to the linear term only 


a 


or 








the water-vapor corrected radiative temperature is therefore given 


accurately by 


el) ece][eu] 


From (15), a first estimate of the surface and/or cloud top 
temperature is obtained based upon the averaged form of the regression 
equation (7). This is justifiable insofar as (13) is concerned since 
the coefficients B, and B> vary from their means by only 5%, an insig- 


nificant deviation from the mean. 
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The bracketed factor in (15) may be written as a depletive 


correction to the interface temperature AT (0) 


AT, (0) = 4. (0) Te Uy 


Pe) 
oN 1] O=0...,60 
rv) 
he (0) eee & ) ’EC 6 
Cain 
Values forQT,,( ©) are listed in Tables 4, 5.......10, for © =0°, 
10°... eee. -. GO, respectively. The input data for each table 


is . CO) and the optical depth Uk. 

Tests of the tables/graphs against known interface tempera- 
tures show very good results. The corrections AT, are less than 
2K, and usually less than 0.5K, but always additive in the sense 
of adding back depleted flux. In fact, the two steps just described, 
namely forming 

Ty + ATs 

results in specification of the interface temperature T accurate to 
+ 0.2K. An alternative method for obtaining A\ Tp is discussed in 


Sectionw- 


7. Angular-Dependence Corrections and Conclusions 

The major results derived have been the determination of a 
composite, pooled, first estimate of 7, from (13). Since i has 
been derived by averaging al! of the nearly identical forms of the 
terms By +B 5X5 of regression equation (7), 1, may be considered 
independent of the zenith angle. 


The depletive correction Al, (2) to the radiative temperature, 


arising from the presence of the water-vapor term was of the form 


ATz (@) = k(8) Te Ue 
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in the regression equation (7, 15). Here k (©) is a constant 
for each selection of O , but differs from one zenith angle to 
the next. 

To emphasize the nature of the angular dependency which now 
resides solely in the depletion term AT, (@) , consider the 
cases O = 00 and 60°: 

@=0° ATB(u,0)=2 S66 Ty lly X/0™ 
= 3 

@=60° AT;(y, SéC60) = 5 2987, BU, X70 
This specialization is done not only because these results span 
completely the angular range employed, but also because both cases 


corresponded to significant F'-tests for variables X. in the analysis 


3 
of variance (Section 4). 
For intermediate values of 8,475 (@) may be estimated by 


the Taylor expansion formula applied in the form 


AT, (@) = OT; ()+[ 27 Jacsec Qe) (16) 





and when finite-differences are employed, 
A7z (6) = ATz (0) + ATs (60)-A7s (2)/(eece-1) 
, 


from which we finally obtain 


“N 
AT; (6)= 4, 7 [2 932sec0- 0.566)x 10 (17) 


(17) fits the results of Appendix II, Tables 4 - 10, within + 0.2K 
and may be considered a good estimator of the best-fit depletive- 


correction for all..angles O in the range of O x 609. 
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Further, the combination of (17) together with 7, of (13) led to 
the recovery of the known interface temperature to an accuracy of 
FuOBZK. 

As a measure of the effects of energy depletion, consider 


the intensity-depletion A Z(6) 


ee ; [ a ) ] “ a hay 


" AT(®) _ 4a7;(@) 


—_ 
= 


Tp ( 7s) a alts) 








With regard to flux absorbed in the atmosphere, relative to the 
m—N 
cloud-top estimate, Fa ( T +£\T), one may integrateA 1¢ ©) 


hemispherically to obtain /p 
A F=/0°F, (7g) un J #2 952 SEC 6-Sbe]sinecosede 
O 


Upon integration, the last equation gives the flux-loss between 


interface and the channel 2 sensor 
e oN ae 
AF = Fy (Tz) 4, [10.57 x/0 >| (19) 


Here the last expression in the bracket shows the result of the 
integration. Limb-darkening has been neglected, since according to 
Wark et al [7] only about 25% of the possible flux lies outside 
of a 60° angle cone. 

Even though the cloud temperature sensed by Channel 2 may be 


too low by 0.5° to 2.0°, depending upon the zenith angle, the 


"sensed" black-body intensity and flux values at cloud tops are too 
low by approximately 1 to 2% depending upon the actual value of 
ATp (0)/Ts of (17) and upon the value Ux in (19), respectively. 
Refinements of this statistical model may be made to include 
a posteriori values of ozone or of dust-scattering, but such 
refinements are not likely to be as statistically significant on 
an operational basis where the optical masses of ozone and of dust 


are not known until long after the data readout is made. 
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APPENDIX I 


Location 
Havana, Cuba 
Kenitra, Morocco 
Oakland, California 
Nantucket, Massachusetts 
Ft. Smith, Arkansas 
Thule, Greenland 
Chaquaromos Bay, Trinidad 
Tucson, Arizona 
Washington, D.C. 
Clark Field, P.I. 
Stephenville, Texas 
Glascow, Montana 
Barter Island, Alaska 
Chateauroux, France 


Bitteurg, Germany 
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1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 
1200 


15 Clear Sky Atmospheric Models from [7 | 


Used in Multiple Linear Regression Program 


GMT , 
GMT , 
GMT , 
GMT , 
GMT , 
GMT , 
GMT , 
GMT, 
GMT , 
GMT , 
GMT , 
GMT , 
GMT, 
GMT , 
GMT, 


Time 


Sept. 
Sept. 
Sept. 
Sept. 
Sept. 
Sept. 


29, 19S 
ey, ee 
29, 1958 
29, 198 
29, 1958 
29, l9e8 


April 1, 1958 


April 1, 1958 


April 1, 1958 


Jan. 


ib, 1968 


April 1, 1958 


Jan. 
Jan. 
Jan. 


Jan. 


1, 1958 
is, 1966 
co, i900 
25 aul 9e 


APPENDIX II 
wa. 
CORRECTIONS TO FIRST ESTIMATES OF Tp BASED 
ON WATER VAPOR PATH FOR ZENITH ANGLESQO=10°, 
32.5°, 459, 52.59, and 60° (Tables 4, 5, 6, 
7, 8, 9 and 10) 
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TABLE 4 


Corrections to First Estimates of ip Based 


aN 0.1 


0 

2 

4 

6 

8 
10 
12 
14 
16 
18 
20 
(Ma 
24 
26 
28 
30 
Se 
34 
36 
38 
40 


.0424 
.0462 
.0493 
,0577 
.0540 
.0561 
mO)5) 7, 
mobepe ls 
.0611 
.0620 
.0638 
.0650 
. 0662 
70673 
.0684 
20605 
.0704 
OE 
iae 
07 Sil 
MOV Ecte 


on Water Vapor Path for Zenith Angle 


OZ 


.0848 
.0924 
.0986 
. 1034 
. 1080 
mle 
. 116s 
- 1190 
222 
. 1240 
. 1276 
. 1300 
. 1324 
. 1346 
. 1368 
1350 
. 1408 
. 1424 
. 1444 
. 1462 
. 1478 


oO 


0.4 


. 1696 
. 1848 
Ve 
. 2068 
.2160 
.2244 
. CSO 
. 2380 
.2444 
. 2480 
peoa2 
. 2600 
. 2648 
.2692 
2136 
. 2780 
.2816 
. 2848 
. 2888 
.2924 
. 2956 


30 


. 3392 
. 3696 
. 3944 
.4126 
. 4320 
. 4488 
- 4632 
-4760 
. 4888 
- 4960 
.5104 
socuU 
.5296 
.5384 
.9472 
.9960 
00sec 
. 9696 
S70 
. 9848 
Beep 


.6784 
sf 382 
. /888 
8272 
. 8640 
. BOR 
.9264 
9520 
BLS 
oa) 
.0208 
.0400 
- 0592 
.0768 
.0944 
» Ve 
. 1264 
. 1392 
PF sieya 
. 1696 
. 1824 


NO RO Hd BRO BRO WR) RO NO NO RO A) ] ] a Ss es aes as es es 


. IGS 
.4784 
S07 TO 
.6544 
. 7280 
. Pee 
.8528 
. 9040 
. 9552 
. 9840 
.0416 
.0800 
. 1184 
. 1536 
. 1888 
. 2240 
.2528 
.3784 
.3104 
. Idee 
. 3648 


TABLE 5 


a 
Corrections to First Estimates of T. Based 


Pe 0.1 


.0467 
.0508 
.0541 
. 0569 
.0594 
. OGM 
.0635 
. 0652 
.0670 
.0685 
.0800 
.O7mis 
.0726 
.0738 
.0750 
.0760 
.0871 
- OFZ 
. 0792 
. 0802 
.081] 


on Water Vapor Path for Zenith An§le 


Oa 


. 0963 
- 1016 
. 1082 
.1138 
. 1188 
. 1268 
. 1270 
. 1306 
. 1340 
. 1370 
. 1400 
. 1426 
. 1452 
. 1476 
. 1500 
. 1520 
. 1542 
. 1564 
. 1584 
. 1604 
. 1622 


6=/0° 


0.4 


. 1866 
. 2082 
.2164 
.22/6 
2070 
2460 
. 2540 
. 26i2 
. 2680 
.2740 
. 2800 
cove 
. 2904 
Coe 
. 3000 
. 3040 
. 3084 
. 326 
.3168 
. 3208 
.3244 
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 Sei2 
. 4064 
4328 
~4552 
-4752 
-4920 
. 9080 
.5224 
- 5360 
. 5480 
. 9600 
.5704 
. 9808 
.5904 
. 6000 
.6080 
.6168 
.6256 
.6336 
6416 
. 6488 


eel eed eed eed od od ed od Oe C8! CO CB C/I CC! C/I 


. 7464 
.8128 
.8656 
-9104 
. 9504 
.9840 
.0160 
.0448 
SOO 
.0960 
. 1200 
. 1408 
. 1616 
. 1808 
. 2000 
.2160 
, 2380 
.2512 
,cOra 
. 2832 
Cone 


NO BRO RS BPO BS POS BS DPS DP PO PO PPD D a A st ss OO SG! 


.4928 
.6256 
. Fam 
.8208 
. 9008 
. 9680 
.0320 
.0896 
. 1440 
. 1920 
. 2400 
.2816 
. 32392 
. SMS 
- 4000 
-4320 
4672 
.5024 
.5344 
. 5664 
9 


TABLE 6 


“N 
Corrections to First Estimates of T, Based 


Uy} 0.1 
L(0 

mo 6' 0314 
2 0346 
. 03% 
=  .039) 
8  .0410 
10 0425 
I .0440 
14 = .0452 
16 .0464 
18  .0476 
20 ~=—«-«.0486 
22 ~=—- 0496 
24 ~=—.0505 
26 ~—S.0513 
Ze | C052) 
30 ~=—-.0530 
32 = .0536 
34.0545 
36 =. 0551 
38. = .0557 
40 .0563 


on Water Vapor Path for Zenith An§le 


OR 


.0628 
. 0692 
.0742 
OFZ 
.0820 
.0850 
.0880 
.0904 
O9ize 
. 09152 
C372 
0992 
. 1010 
. 02s 
. 1042 
. 1060 
OWE 
. 1090 
.1102 
.1114 
.1126 


9=20° 


0.4 


. 1256 
. 1384 
. 1484 
. 1564 
. 1640 
. 1700 
. 1760 
. 1808 
. 1856 
. 1904 
. 1944 
. 1984 
. 2020 
. 20872 
. 2084 
. 2020 
.2144 
.2180 
.2204 
cea 
Say 


ae 


Ome 


262 
.2/768 
. 2968 
.3128 
. 3280 
. 3400 
. 3920 
. 3616 
3) We 
. 3808 
. 3888 
. ICS 
. 4040 
.4104 
.4168 
.4240 
- 4288 
. 4360 
-4408 
-4456 
.4504 


.5024 
. 5536 
. 5986 
.6256 
.6560 
. 6800 
. 7040 
. 7262 
.7424 
. 7616 
7196 
. 1936 
. 8080 
.8208 
. 8386 
. 8480 
.8576 
.8720 
.8826 
.8912 
. 9008 


SZ 


.0048 
. OFZ 
. 1872 
. 2en2 
.3120 
. 3600 
- 4080 
-4464 
. 4848 
. amie 
. 5952 
. 5OFZ 
-4160 
-6416 
. 66m 
.6960 
. GZ 
. 7440 
. (OZ 
. 7824 
.8016 


TABLE 7 


oN 
Corrections to First Estimates of Ts Based 
on Water Vapor Path for Zenith AnBle 


O=Dhug 
Te) a 0.1 0.2 0.4 0.8 1.6 22 

0 .0508 .1016 . 2032 .4064 .8128 1.6256 
me =. .0550 .1100 .2200 .4400 . 8800 1.7600 
4 0586 alg .2344 . 4688 . 9376 1.8752 
om 10615 .1230 . 2460 .4920 . 9840 1.9680 
8 = .0641 . 1282 . 2564 sZeR | 50256 2 ORTIZ 
10 = .0663 . 1326 .2652 5304 = 1.0608 2 . VaiNG 
12.0685 . 1370 . 2740 5480 1.0960 2.1920 
14.0704 . 1408 . 2816 5632 1.1264 2.2528 
ver 0722 1444 . 2888 6776" ‘laleee 2.3104 
fee = 0737 1474 . 2948 5896 Val 92 2.3584 
20 = .0753 . 1506 .3012 6024 1.2048 2.4096 
a OnGe . 1536 .3072 6144 1.2288 2.4576 
24.0781 . 1562 .3124 6248 1.2496 2.4992 
Ze OVS . 1586 rol/2 6344 =: 1.2688 2.5376 
28  — .0807 .1614 .3228 6456 1.2912 2.5824 
30 = .0818 . 1636 272 6544 1.3088 = 2.6276 
32s «0829 . 1658 .3316 6632 1.3264 2.6528 
34.0841 . 1682 . 3366 6732 1.3464 2.6928 
36 = .0852 .1704 . 3408 6816 1.3632 2.7264 
38 = .0862 .1724 .3448 6896 1.3792 2.7584 
40 —.0872 .1744 . 3488 “6o76  Ireoo? 2.7904 


so). 


TABLE 8 


oN 
Corrections to First Estimates of T, Based 


T(e) Us| 0.1 
0 .0584 
2 .0639 
4 .0683 
6 0726 
8 0752 

10 .0780 
2 .0806 
14 .0828 
16 .0849 
18 . 0869 
20 .0888 
22 .0906 
24 .0923 
26 . 0938 
28 .0ga5 
30 . C867 
32 0282 
34 .0994 
36 . 1007 
38 1S 
40 . 1031 


on Water Vapor Path for Zenith Angle 


. 1168 
. 1275 
. 1366 
. 1440 
. 1504 
1900 
21612 
J 1655 
. 1698 
ai icte 
mera) 
. Wii 
. 1846 
. 1876 
. 1906 
. 1934 
. 1964 
oe 
.2014 
208s 
. 2062 


0-45° 


0.4 


. 2336 
:2000 
2 PSe 
. 2880 
. 3008 
. 3220 
3229 
3512 
. 336 
. 3476 
3652 
. 3624 
S072 
. 3/52 
eel 2 
. 3868 
. 92S 
. Same 
. 4028 
.4076 
-4124 


0.8 


.4672 
Smiliz 
.5464 
. 5760 
.6016 
.6240 
.6448 
.6624 
OV o2 
6852 
.7104 
. 7248 
. 7384 
. 7504 
. 7624 
. L3G 
a isholo 
. (982 
. 8056 
.8152 
.8248 


.9344 
.0224 
.0928 
1520 
05a 
. 2480 
. 2696 
. 3248 
. 3584 
. 3904 
. 4208 
- 4496 
.4768 
. 5008 
.5248 
5472 
d/h 
.5904 
6112 
.6304 
.6496 


WW WW WW WW PDP MH POP MPP PPP DY — 


CIR 


. 8688 
.0448 
. 1856 
. 3040 
. 4064 
. 4960 
. 582 
.6496 
. 71168 
. 7808 
.8416 
. 8982 
. Fao 
.0016 
.0496 
.0944 
. 1424 
. 1808 
2224 
. 2608 
CFL 


TABLE 9 


Vs 
Corrections to First Estimates of Tp Based 


on Water Vapor Path for Zenith Angle 


0777 
.0852 
.0909 
.0956 
.0999 
. 1034 
. 1068 
. 1100 
.1128 
a5 
myo 
sneO | 
mee 
. 2343 
. 1263 
alco) 
. 1301 
Se 
30 
1390 
. 1366 


ORZ 


. 1554 
. 1704 
. 1818 
leila 
. 1998 
. 2068 
ras 
.2200 
Cale 
52910 
F2308 
.2402 
- 2446 
. 2486 
NeDcO 
2a e¥2 
Ae 
. 2634 
. 26/0 
.2/00 
BOL SC 


B=o2o) 
0.4 0.8 
73106 .6216 
.3408 .6816 
3636 RMF Ma 
3824 ./668 
. 3996 .7992 
.4136 ROZie 
4272 ~8544 
~4400 . 8800 
~4512 ~9024 
~4620 ~9220 
.4716 9432 
~4804 .9608 
~4992 .9874 
~4972 .9944 
~5052 1.0104 
.5124 1.0248 
~5204 1.0408 
.5268 1.0536 
~5340 1.0680 
5400 1.0800 
5464 1.0928 


IPs 


1.2432 
3632 
1.4544 
eoss6 
1.5884 
1.6544 
1.7088 
1. 7600 
1.8048 
1.8440 
1.8864 
A SIeie 
egoGe 
1.9888 
2.0208 
2.0496 
2.0816 
Araley | 
Ze soU 
2.1600 
2.1856 


PPP PP HP BPW WW WW WW WW W W PO PP PO 


SZ 


. 4864 
. 7264 
. 9088 
.0672 
. 1768 
. 3088 
~4176 
.0c0U 
.6096 
. 6880 
11@6 
8432 
EwiSie 
a Tile 
.0416 
0992 
. 1632 
2144 
2720 
. 3200 
Beye || 2 


TABLE 10 


“N 
Corrections to First Estimates of Tp Based 
on Water Vapor Path for Zenith Angle 


oO 
8 =60 
TON 0.7 0.2 0.4 0.8 16 3.2 

0  .0989 1979 3957 7914 1.5828 3.1657 
2 1072 2145 4290 8580 1.7159 Me eas 
Oe ice 2278 4555 9111 1.8222 3.6444 
@ 1192 2395 4791 9582 1.9164 3.8327 
Se  WaaG 2492 4983 9966 1.9932 3.9865 
10.1292 2584 ‘5169 " 1m0@388 ©) 2. leveu. INGEN 
la. .teee 2662 5327 1ap653 2ite06 My. 2082 
ie Ma SGe 2738 5476 «1.0951 2.1902 4.3805 
16  .1404 2808 5616 1.1232 2.2464 4.4928 
18.1434 2867 5735 1.1470 2.2039 @.5a7¢ 
20». 1463 2927 5854 1.1707 = 2.9414 @.6eas 
22 =. 1494 2988 5977 981953 293007 aR 7and 
94 ~—s«. 1518 3037 6074 1.2148 2.4296 4.9591 
26s 1547 3093 6186 1.2372 2.4745 4.9490 
28 =, 1568 Fes 273 «1.2545 |) 2-500) ae MON OMeT 
30 =. 1594 3139 6374 1.2748 2.5497 5.0se6 
NGG 3231 (0463 1.2925 2.5851 6.17@2 
Zim 1635 3269 nepe8 1.0077 296153 Me) 2u0n 
Some | le57 eS 6629 1.3258 2.6516 5.3032 
38—~C«#w« 67 2853 6707 +1.3414 ~=2.6827. «5.3654 
40. 1684 3368 6737. -:1.3474 += 2.6948 ~=—s«455 3896 
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